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Abstract 

The centrality dependence of the mean charged-particle multiplicity as a function of pseu¬ 
dorapidity is measured in approximately 1 pb _1 of proton-lead collisions at a nucleon- 
nucleon centre-of-mass energy of = 5.02 TeV using the ATLAS detector at the Large 
Hadron Collider. Charged particles with absolute pseudorapidity less than 2.7 are recon¬ 
structed using the ATLAS pixel detector. The p+Pb collision centrality is characterised by 
the total transverse energy measured in the Pb-going direction of the forward calorimeter. 
The charged-particle pseudorapidity distributions are found to vary strongly with central¬ 
ity, with an increasing asymmetry between the proton-going and Pb-going directions as the 
collisions become more central. Three different estimations of the number of nucleons parti¬ 
cipating in the p+Pb collision have been earned out using the Glauber model as well as two 
Glauber-Gribov inspired extensions to the Glauber model. Charged-particle multiplicities 
per participant pair are found to vary differently for these three models, highlighting the im¬ 
portance of including colour fluctuations in nucleon-nucleon collisions in the modelling of 
the initial state of p+Pb collisions. 
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1 Introduction 

Proton-nucleus (p+ A) collisions at the Large Hadron Collider (LHC) [1] provide an opportunity to probe 
the physics of the initial state of ultra-relativistic heavy-ion (A+A) collisions without the presence of 
thermalisation and collective evolution [2]. In particular, p+A measurements can provide insight into the 
effect of an extended nuclear target on the dynamics of soft and hard scattering processes and subsequent 
particle production. Historically, measurements of the average charged-particle multiplicity as a function 
of pseudorapidity, dA^h/dr/. where pseudorapidity is defined as tj = - In tan(d/2) with 6 the particle angle 
with respect to the beam direction, have yielded important insight into soft particle production dynamics in 
proton- and deuteron-nucleus ( p/d+A ) collisions [3-8] and provided essential tests of models of inclusive 
soft hadron production. 

Additional information is obtained if measurements of the charged-particle multiplicities are presented 
as a function of centrality, an experimental quantity that characterises the p/d+A collision geometry. 
Previous measurements in d+Au collisions at the Relativistic Heavy Ion Collider (RHIC) [9] have char¬ 
acterised the centrality using particle multiplicities at large pseudorapidity, either symmetric around mid¬ 
rapidity [10] or in the Au fragmentation direction [11]. These measurements have shown that the rapidity- 
integrated particle multiplicity in d +Au collisions scales with the number of inelastically interacting, or 
“participating", nucleons, A part . This scaling behaviour has been interpreted as the result of coherent 
multiple soft interactions of the projectile nucleon in the target nucleus, and is known as the wounded- 
nucleon (WN) model [12]. The charged-particle multiplicity distributions as a function of pseudorapidity 
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measured in central r/+Au collisions are asymmetric and peaked in the Au-going direction [7]. This 
observation has been explained using well-known phenomenology of soft hadron production [13]. 

There are alternative descriptions of the centrality dependence of the d/V c n/d// distribution in r/+Au col¬ 
lisions at RHIC [14, 15] and p+ Pb collisions at the LHC [15-17] based on parton saturation models. 
Measurements of the centrality dependence of dA t h/d// distributions in /;+Pb collisions provide an essen¬ 
tial test of soft hadron production mechanisms at the LHC. Such tests have become of greater importance 
given the observation of two-particle [18-21] and multi-particle [21-23] correlations in the final state of 
p+ Pb collisions at the LHC. These correlations are currently interpreted as resulting from either initial- 
state saturation effects [15, 24, 25] or from the collective dynamics of the final state [26-30]. For either 
interpretation, information on the centrality dependence of d/V c h /d// can provide important input for de¬ 
termining the mechanism responsible for these structures. 

Recent measurements from the ALICE experiment [31] show behaviour in the centrality dependence of 
the charged-particle pseudorapidity distributions, which is qualitatively similar to that observed at RHIC. 
That analysis compared different methods for characterising centrality and suggested that the method used 
to define centrality may have a significant impact on the centrality dependence of the measured dN c h/dp 
distribution. 

An important component of any centrality-dependent analysis is the geometric model used to relate ex¬ 
perimental observables to the geometry of the nuclear collision. Glauber Monte Carlo models [32], which 
simulate the interactions of the incident nucleons using a semi-classical eikonal approximation, have been 
successfully applied to many different A+A measurements at RHIC and the LHC. A key parameter of 
such models is the inelastic nucleon-nucleon cross-section, which is taken to be 70 mb for this ana¬ 
lysis [31]. However, the Glauber multiple-scattering approximation assumes that the nucleons remain on 
the mass shell between successive scatterings, and this assumption is badly broken in ultra-relativistic 
collisions. Corrections to the Glauber model [33], hereafter referred to as “Glauber-Gribov,” are needed 
to account for the off-shell propagation of the nucleons between collisions. 

A particular implementation of the Glauber-Gribov approach is provided by the colour-fluctuation model 
[34-37], That model accounts for event-to-event fluctuations in the configuration of the incoming proton 
that are assumed to be frozen over the timescale of a collision and that can change the effective cross- 
section with which the proton scatters off nucleons in the nucleus. These event-by-event fluctuations in the 
cross-section can be represented by a probability distribution P(cr). The width of that distribution can be 
characterised by a parameter uv, which is the relative variance of the cr distribution, u> (T = ((cr/cr to t - l) 2 ). 
The usual total cross-section, cr tot , is the event-averaged cross-section, or, equivalently, the first moment of 
the P(cr) distribution, cr tot = J () dcrP(cr ) cr. The parameter oj <r can be measured using diffractive proton- 
proton scattering at high energy [35, 36]. First estimates of ov at LHC energies [36] extrapolated to 5 TeV 
yielded oJa- ~ 0.11, while a more recent analysis suggested ov ~ 0.2 [37]. Because the cross-section 
fluctuations in the Glauber-Gribov colour-fluctuation (GGCF) model may have a significant impact on 
the interpretation of the results of this analysis, the geometry of p+ Pb collisions has been evaluated using 
both the standard Glauber model and the GGCF model with to tT = 0.11 and 0.2. 

This paper presents measurements of the centrality dependence of dN C h/dp in p+ Pb collisions at = 
5.02 TeV using 1 pb -1 of data recorded by the ATLAS experiment [38] in September 2012. Charged 
particles are detected in the ATLAS pixel detector and are reconstructed using a two-point tracklet al¬ 
gorithm similar to that used for the Pb+Pb multiplicity measurement [39]. Measurements of dN^/dp are 
presented for several intervals in collision centrality characterised by the total transverse energy measured 
in the forward section of the ATLAS calorimeter on the Pb-going side of the detector. A standard Glauber 
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model [32] and the GGCF model [36, 37] with oj (T = 0.11 and 0.2 arc used to estimate (V part ) for each 
centrality interval, allowing a measurement of the /V part dependence of the charged-particle multiplicity. 

The paper is organised as follows. Section 2 describes the subdetectors of the ATLAS experiment rel¬ 
evant for this measurement. Section 3 describes the event selection. Section 4 describes the Monte 
Carlo simulations used to understand the performance and derive the corrections to the measured quant¬ 
ities. Section 5 describes the choice of centrality variable. Section 6 describes the measurement of the 
charged-particle multiplicity and Sect. 7 describes the estimation of the systematic uncertainties. Sec¬ 
tion 8 presents the results of the measurement, and the interpretation of the yields of charged particles 
per participant is discussed in Sect. 9. Section 10 concludes the paper. The estimation of the geomet¬ 
ric parameters in each centrality interval for the Glauber and GGCF models is presented in detail in the 
Appendix. 


2 Experimental setup 


The ATLAS detector is described in detail in Ref. [38]. The data selection and analysis presented in this 
paper is performed using the ATLAS inner detector (ID), calorimeters, minimum-bias trigger scintillators 
(MBTS), and the trigger system. The inner detector measures charged-particle tracks using a combination 
of silicon pixel detectors, silicon microstrip detectors (SCT), and a straw-tube transition-radiation tracker 
(TRT), all immersed in a 2 T axial magnetic field. The pixel detector is divided into “barrel” and “endcap” 
sections. For collisions occurring at the nominal interaction point, 1 the barrel section of the pixel detector 
allows measurements of charged-particle tracks over |//| < 2.2. The endcap sections extend the detector 
coverage, spanning the pseudorapidity interval 1.6 < |//| < 2.7. The SCT and TRT detectors cover 
|//| < 2.5 and |//| < 2, respectively, also through a combination of barrel and endcap sections. 

The barrel section of the pixel detector consists of three layers of staves at radii of 50.5 mm, 88.5 mm, 
and 122.5 mm from the nominal beam axis, and extending ±400.5 mm from the centre of the detector in 
the z direction. The endcap consists of three disks placed symmetrically on each side of the interaction 
region at z locations of ±493 mm, ±578 mm and ±648 mm from the centre of the detector. All pixel 
sensors in the pixel detector, in both the barrel and endcap regions, are identical and have a nominal size 
of 50/im x 400 pm. 

The MBTS detect charged particles in the range 2.1 < |//| < 3.9 using two hodoscopes, each of which is 
subdivided into 16 counters positioned at z = ±3.6 m. The ATLAS calorimeters cover the full azimuth 
and the pseudorapidity range |//| < 4.9 with the forward part (FCal) consisting of two modules positioned 
on either side of the interaction region and covering 3.1 < |//| < 4.9. The FCal modules are composed 
of tungsten and copper absorbers with liquid argon as the active medium, which together provide 10 
interaction lengths of material. 

The LF1C delivered its first proton-nucleus collisions in a short /;+Pb “pilot” run at y i v NN = 5.02 TeV in 
September 2012. During that run the LHC was configured with a clockwise 4 TeV proton beam and an 
anti-clockwise 1.57 TeV per-nucleon Pb beam that together produced collisions with a nucleon-nucleon 
centre-of-mass energy of ^Js SN = 5.02 TeV and a longitudinal rapidity boost of 0.465 units with respect to 


1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector 
and the z-axis along the beam pipe. The .r-axis points from the IP to the centre of the LHC ring, and the (/-axis points 
upward. Cylindrical coordinates (r, <f>) are used in the transverse plane, 0 being the azimuthal angle around the beam pipe. 
The pseudorapidity is defined in terms of the polar angle 6 as rj = - In tan(0/2). 
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the ATLAS laboratory frame. Following a common convention used for p+A measurements, the rapidity 
is taken to be positive in the direction of the proton beam, i.e. opposite to the usual ATLAS convention for 
pp collisions. With this convention, the ATLAS laboratory frame rapidity y and the p+Pb centre-of-mass 
system rapidity y cm are related as y cm = y - 0.465. 


3 Event selection 

Minimum-bias p+Pb collisions were selected by a trigger that required a signal in at least two MBTS 
counters. The p+Pb events selected for analysis are required to have at least one hit in each side of 
the MBTS, a difference between the times measured in the two MBTS hodoscopes of less than 10 ns, 
and a reconstructed collision vertex in longitudinal direction, z vtx , within 175 mm of the nominal centre 
of the ATLAS detector. Collision vertices are defined using charged-particle tracks reconstructed by 
an algorithm optimised for pp minimum-bias measurements [40]. Reconstructed vertices are required 
to have at least two tracks with transverse momentum pj > 0.4 GeV. Events containing multiple p+Pb 
collisions are rare due to very low instantaneous luminosity during the pilot run and are further suppressed 
in the analysis by rejecting events with two collision vertices that are separated in z by more than 15 mm. 
Applying this selection reduces the fraction of events with multiple collisions from less than 0.07% to 
below 0.01%. 

To remove potentially significant contributions from electromagnetic and diffractive processes, the topo¬ 
logy of the events was first analysed in a manner similar to that performed in a measurement of rapidity 
gap cross-sections in 7 TeV proton-proton collisions [41]. The pseudorapidity coverage of the calori¬ 
meter, -4.9 < j] < 4.9, is divided into A;/ = 0.2 intervals, and each interval containing one or more 
clusters with px greater than 0.2 GeV is considered as occupied. To suppress the contributions from 
noise, clusters are considered only if they contained at least one cell with an energy at least four times the 
standard deviation of the cell noise distribution. 

Then, the edge-gap on the Pb-going side of the detector is calculated as the distance in pseudorapidity 
between the detector edge rj = -4.9 and the nearest occupied interval. Events with edge-gaps larger than 
two units of pseudorapidity typically result from electromagnetic or diffractive excitation of the proton 
and are removed from the analysis. The effect of this selection is identical to the requirement of a cluster 
with transverse energy Ej > 0.2 GeV to be present in the region // < -2.9. No requirement is imposed 
on edge-gaps on the proton-going side. The gap requirement removes, with good efficiency, a sample 
of events which are not naturally described in a Glauber picture of p+Pb collisions. This requirement 
removes a further fraction / gap = 1% of the events passing the vertex and MBTS cuts, yielding a total 
of 2.1 million events for this analysis. The result of this event selection is to isolate a fiducial class of 
p+Pb events, defined as inelastic p+Pb events that have a suppressed contribution from diffractive proton 
excitation events. 

4 Monte Carlo simulation 

The response of the ATLAS detector and the performance of the charged-particle reconstruction al¬ 
gorithms are evaluated using one million minimum-bias 5.02 TeV Monte Carlo (MC) p+Pb events, pro¬ 
duced by version 1.38b of the Hijing event generator [42] with diffractive processes disabled. The four- 
momentum of each generated particle is longitudinally boosted by a rapidity of 0.465 to match the beam 
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conditions in the data. The detector response to these events is fully simulated using Geant4 [43, 44]. The 
resulting events are digitised using conditions appropriate for the pilot /j+Pb run and fully reconstructed 
using the same algorithms that are applied to the experimental data. This MC sample is primarily used to 
evaluate the efficiency of the ATLAS detector for the charged-particle measurements. 

The detector response and event selection efficiencies for peripheral and diffractive p +Pb events have 
properties similar to those for inelastic or diffractive pp collisions, respectively. To evaluate these re¬ 
sponses and efficiencies, the pp samples are generated at sfs = 5.02 TeV with particle kinematics boosted 
to match the /;+Pb beam conditions. Separate samples of minimum-bias, single-diffractive, and double- 
diffractive pp collisions with one million events each are produced using both Pythia6 [45] (version 
6.425, AMBT2 parameter set (tune) [46], CTEQ6L1 PDF [47]) and Pythia8 [48] (version 8.150, 4C 
tune [49], MSTW2008LO PDF [50]), and simulated, digitised and reconstructed in the same manner as 
the p +Pb events. These six samples are primarily used for the Glauber model analysis described in the 
Appendix. 


5 Centrality selection 


For Pb+Pb collisions, the ATLAS experiment uses the total transverse energy, measured in the two 
forward calorimeter sections to characterise the collision centrality [51]. However, the intrinsic asym¬ 
metry of the p +Pb collisions and the rapidity shift of the centre-of-mass causes an asymmetry in the 
soft particle production measured on the two sides of the calorimeter. Figure 1 shows the correlation 
between the summed transverse energies measured in the proton-going (3.1 < rj < 4.9) and Pb-going 
(-4.9 < T] < -3.1) directions, X £■{- and XA l | ' b , respectively. The transverse energies are evaluated at an 
energy scale calibrated for electromagnetic showers and have not been corrected for hadronic response. 
Figure 1 shows that the mean X Ej rapidly flattens with increasing Xfor > 30 GeV, indic- 
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Figure 1: Distribution of proton-going (YjEj) versus Pb-going (E-Eij?) total transverse energy in the forward calo¬ 
rimeter for p+Pb collisions included in this analysis. The curve shows the average £ E!’ as a function of 2 /i l | >h . 

ating that YjEj i s l ess sensitive than YjEj b to the increased particle production expected to result from 
multiple interactions of the proton in the target nucleus in central collisions. Thus, X Ej alone, rather 
than Xis chosen as the primary quantity used to characterise p+Pb collision centrality for the 
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measurement presented in this paper. However, we describe alternate choices of the centrality-defining 
region below and evaluate the sensitivity of the measurement to this definition. 

The distribution of * ()l events passing the p+Pb analysis selection is shown in Fig. 2. The fol¬ 
lowing centrality intervals are defined in terms of percentiles of the distribution: 0-1%, 1-5%, 

5-10%, 10-20%, 20-30%, 30-40%, 40-60%, and 60-90%. The 2fs!p b ran S cs corresponding to these 
centrality intervals are indicated by the alternating filled and unfilled regions in Fig. 2, with the 0-1% 
interval, containing the most central collisions, being rightmost. Since the composition of the events in 
the most peripheral 90-100% interval is not well constrained, these events are excluded from the ana¬ 
lysis. The nominal centrality intervals were defined after accounting for a 2% inefficiency, as described 
in the Appendix, for the fiducial class of p+Pb events defined above to pass the applied event selection. 
Alternate intervals were also defined by varying this estimated inefficiency to 0% and 4%, and is used 
as a systematic check on the results. While the inefficiency is confined to the 90-100% interval, it in¬ 
fluences the ranges associated with each centrality interval. Potential hard scattering contributions 
to X f2|- b have been evaluated in a separate analysis [52] by explicitly subtracting the contributions from 
reconstructed jets that fall partly or completely in the Pb-going FCal acceptance. That analysis showed 
negligible impact from hard scattering processes on the measured 2^x b distribution. 



Figure 2: Distribution of the Pb-going total transverse energy in the forward calorimeter values for events 

satisfying all analysis cuts including the Pb-going rapidity gap exclusion. The alternating shaded and unshaded 
bands indicate centrality intervals, from right (central) to left (peripheral), 0-1%, 1-5%, 5-10%, 10-20%, 20-30%, 
30-40%, 40-60%, 60-90% and the interval 90-100% that is not used in this analysis. 

To test the sensitivity of the results to the choice of pseudorapidity interval used for the measurement, 
two alternative quantities are defined. The former, £ E ^ , is defined as the total transverse energy 
in FCal cells with // < -4.0. The latter, j s defined as the total transverse energy in the 

two intervals 4.0 < 77 < 4.9 and -4.0 < // < -3.1, an approximately symmetric interval when expressed 
in pseudorapidity in the centre of mass system rj cm . The first of these alternatives is used to evaluate 
the potential auto-correlation between the measured charged-particle multiplicities and the centrality ob¬ 
servable by increasing the rapidity gap between the two measurements. The second is used to evaluate 
the differences between an asymmetric (Pb-going) and symmetric (both sides) centrality observable. The 
effect of these alternative definitions is discussed in Sect. 8. 

The Glauber analysis [32] was applied to estimate (Apart) for each of the centrality intervals used in this 
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analysis. A detailed description is given in the Appendix; only a brief summary of the method is given 
here. The PHOBOS MC program [53] was used to simulate the geometry of inelastic p +Pb collisions 
using both the standard Glauber and GGCF models. The resulting /V part distributions are convolved with 
a model of the /V part -dependent £fd! b distributions, the parameters of which are obtained by fitting the 
measured J]Ej b distribution. The average A part associated with each centrality interval is obtained with 
systematic uncertainties. The results are shown in Fig. 3 for the Glauber model and for the GGCF model 
with (jJq- = 0.11 and 0.2. 
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Figure 3: Mean value of the number of participating nucleons (A part ) for different centrality bins, resulting from fits 
to the measured distribution using Glauber and Glauber-Gribov A part distributions. The error bars indicate 

asymmetric systematic uncertainties. 


6 Measurement of charged-particle multiplicity 

6.1 Two-point tracklet and pixel track methods 

The measurement of the charged-particle multiplicity is performed using only the pixel detector to max¬ 
imise the efficiency for reconstructing charged particles with low transverse momenta. Two approaches 
are used in this analysis. The first is the two-point tracklet method commonly used in heavy-ion collision 
experiments [39, 54, 55]. Two variants of this method are implemented in this analysis to construct the 
d/Vch/d T] distribution and to estimate the systematic uncertainties, as described below. The second method 
uses “pixel tracks”, obtained by applying the full track reconstruction algorithm [56] only to the pixel de¬ 
tector. The pixel tracking is less efficient than the tracklet method as is justified later in the text, but 
provides measurements of the particle pj. The dN c h/dij distribution measured using pixel tracks provides 
a cross-check on the primary measurement that is performed using the two-point tracklets. 

In the two-point tracklet algorithm, the event vertex and clusters [57] on an inner pixel layer define a 
search region for clusters in the outer layers. The algorithm uses all clusters, except the clusters which 
have low energy deposits inconsistent with minimum-ionising particles originating from the primary ver¬ 
tex. The algorithm also rejects duplicate clusters resulting from the overlap of the pixel sensors or arising 




from a small set of pixels at the centre of the pixel modules that share readout channels [58]. Two 
clusters in a given layer of the pixel detector are considered as one if they have an angular separation 
VwF" + (dp) 2 < 0 . 02 , because they likely result from the passage of a single particle. 

The pseudorapidity and azimuthal angle of the cluster in the innermost layer (p, (p) and their differences 
between the outer and inner layers (A p, A <p) are taken as the parameters of the reconstructed tracklet. The 
A p of a tracklet is largely determined by the multiple scattering of the incident particles in the material 
of the beam pipe and detector. This effect plays a less significant role in the A (p of a tracklet, which is 
driven primarily by the bending of charged particles in the magnetic field, and hence one expects A (p to 
be larger. The tracklet selection cuts are: 

|A 17 I < 0.015, \A(p\ < 0.1, (1) 

|Ap| < \A(p\. ( 2 ) 

Keeping tracklets with \ A(p\ <0.1 corresponds to accepting particles with pj >0.1 GeV. The selection in 
Eq. (2) accounts for the momentum dependence of charged-particle multiple scattering. 

The Monte Carlo simulation for the dA c h /d /7 analysis is based on the Hijing event generator, which is 
described in Sect. 4. The Hijing event generator is known to not accurately reproduce the measured 
particle pj distributions. This is addressed by reweighting the Hijing p\ distribution using the ratio 
of reconstructed spectra measured with the pixel track method in the data and in the MC simulation. 
The reweighting function is extrapolated below pj = 0.1 GeV and applied to all generated particles 
and their decay products. This is done in intervals of centrality and pseudorapidity. Generator-level 
primary particles are defined as particles with a mean lifetime r > 0.3 • 10 1(1 s either directly produced 
in p+ Pb interactions or from subsequent decays of particles with a shorter lifetime. This definition is 
the same as used in previous measurements of charged-particle production in pp [40] and Pb+Pb [59] 
collisions by ATLAS. All other particles are defined as secondaries. Tracklets are classified as primary or 
secondary depending on whether the associated generator-level charged particle is primary or secondary. 
Association between the tracklets and the generator-level particles is based on the Geant4 information 
about hits produced by these particles. Tracklets that are formed from the random association of hits 
produced by unrelated particles, or hits in the detector which are not matched to any generated particle 
are referred to as “fake” tracklets. 

The contribution of fake tracklets is relatively difficult to model in the simulation, because of the a priori 
unknown contributions of multiple sources, such as noisy clusters or very low energy particles. To address 
this problem, the tracklet algorithm is used in two different implementations referred to as “Method 1” 
and “Method 2”. In Method 1, at most one tracklet is reconstructed for each cluster on the first pixel 
layer. If multiple clusters on the second pixel layer fall within the search region, the resulting tracklets 
are merged into a single tracklet. This approach reduces, but does not eliminate, the contribution of fake 
tracklets that are then accounted for using an MC-based correction. Method 2 reconstructs tracklets for 
all combinations of clusters in only two pixel layers, the innermost and the next-to-innermost detector 
layers. To account for the fake tracklets arising from random combinations of clusters, the same analysis 
is performed after inverting the x and y positions of all clusters on the second layer with respect to the 
primary vertex (x - x wtK , y - y vtx ) —> (-(x - x vtx ), -(y - y vtx )). The tracklet yield from this “flipped” 
analysis, A®, is then subtracted from the original tracklet yield, A t e r v to obtain an estimated yield of true 
tracklets N u , 

tftrto) = Kiv) - Kip)- (3) 
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Distributions of A r\ and Ad> of reconstructed tracklets using Method 1 for data and simulated events arc 
shown in Fig. 4 for the barrel (upper plots) and endcap (lower plots) parts of the pixel detector. The 



Figure 4: Stacked histograms for the differences between the hits of the tracklet in outer and inner detector layers 
in pseudorapidity A?/ (left) and in azimuth At/> (right) for the tracklets reconstructed with Method 1 measured in 
the data (points) and simulation (histograms) in p+Pb collisions at -\[s^ =5.02 TeV for barrel (top) and endcaps 
(bottom). Contributions from primary, secondary, and fake tracklets in the simulation are shown separately. The 
lower panels show the ratio of the simulation to the data. 


simulation results show the three contributions from primary, secondary and fake tracklets. The selection 
criteria specified by Eq. (1) are shown in Fig. 4 as vertical lines and applied in A (p for A;/ plots and vice 
versa. Outside those lines, the contributions from secondary and fake tracklets are more difficult to take 
into account, especially in the endcap region. These contributions partially arise from low-px particles on 
spiral trajectories and their description in the MC simulation is therefore very sensitive to the amount of 
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detector material. The ratio between simulation and the data is also shown for each plot. These ratios are 
closer to unity in the barrel region than in the endcap region, where they deviate by up to 5% except at 
very low |A0|. At low |A0| corresponding to high pj, the MC deviates from the data even after reweighing 
procedure based on pixel tracks. This is due to low resolution of pixel track at high pj, however, the 
contribution of high-px particles to dN c h/dij is negligible. 

The top left panel of Fig. 5 shows the pseudorapidity distribution of tracklets reconstructed with Method 2 
and satisfying the criteria of Eqs. (1) and (2) in the 0-10% centrality interval for data (markers) and for the 
simulation (lines). The results of flipped reconstruction are also shown in the plot. The direct and flipped 
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Figure 5: ^-distribution of the number of tracklets reconstructed with Method2. Left top panel: comparison of 
the simulation (lines) to the data (markers). The results of the flipped reconsttuction are shown with open markers 
for data and dashed line for simulation. Right top panel: the simulated result for three contributions: primary, 
secondary and fake tracklets. Square markers show the result of simulation obtained with flipped reconstruction 
events. Lower panels: on the left are the ratios of flipped (A®) to direct (A^ v ) distribution in the data (markers) and 
in the simulation (dashed line); on the right is the ratio of the number of fake tracklets to the number of flipped 
tracklets. 


distributions are each similar between data and MC simulation but not identical, reflecting the fact that 
Hijing does not reproduce the data in detail. Flowever, the lower panel of Fig. 5 shows that the ratios of the 
distribution of the number of tracklets in the flipped and direct events are very similar between the data 
and the MC simulation. A breakdown of the MC simulation distribution into primary, secondary and fake 
tracklets contributions is shown in the top right panel of Fig. 5. The distribution of plotted with 

open markers, closely follows the histogram of fake tracklets. The lower panel of the plot shows the ratio 
of fake tracklet distribution to the flipped distribution. This ratio is consistent with unity to within 5% in 
the entire range of measured ?/. This agreement justifies the subtraction of the fake tracklet contribution 
according to Eq. (3) for Method 2. 

Although the fake rate is largest in Method 2, the flipped method is used to estimate the rate directly 
from the data. In the 0-10% centrality interval, the fake tracklet contribution estimated with this method 
amounts to 8% of the yield at mid-pseudorapidity and up to 16% at large pseudorapidity. In the same 
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centrality interval, the fake tracklet contributions using Method 1 and the pixel track method are smaller, 
vary from 2% to 10% and 0.2% to 1.5%, respectively, but are determined with MC. All three methods 
rely on the MC simulation to correct for the contribution of secondary particles. 
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Figure 6: Left top: distribution from the MC simulation for the generated number of primary charged particles per 
event (tl/V dl /d//) shown with line, reconstructed number per event (1 /N evt dN n -/dri) of track lets from Method 1 shown 
with circles, tracklets from Method 2 after flipped event subtraction shown with squares, and pixel tracks shown 
with diamonds. Left bottom: the ratio of reconstructed to generated tracklets and pixel tracks. Right top: open 
markers represent the same 1/7V ev tdiV tl /d;/ distributions as in the left panel, reconstructed in the data. Filled markers 
of the same shape represent corrected distributions corresponding to diV c h/d//. Right bottom: The ratio of corrected 
distributions of Method 2 and pixel tracks to Method 1. 


6.2 Extraction of the charged-particle distribution 


The data analysis and corresponding corrections are performed in eight intervals of detector occupancy 
(O) parameterised using the number of reconstructed clusters in the first pixel layer and chosen to cor¬ 
respond to the eight p+Pb centrality intervals, and in seven intervals of z vtx , each 50 mm wide. For each 
analysis method, a set of multiplicative correction factors is obtained from MC simulations according to 


C(<9,z vtx ,//) 


Apr(f2, Zvtx 1 1]) 
Arec(f5, Zvtx> V) 


(4) 


Here, N pr and /V rec represent the number of primary charged particles at the generator level and the number 
of tracks or tracklets at the reconstruction level, respectively. These correction factors account for several 
effects: inactive areas in the detector and reconstruction efficiency, contributions of residual fake and 
secondary particles, and losses due to track or tracklet selection cuts including particles with pj below 
0.1 GeV. They are evaluated as a function of (9, z v tx, and rj both for the fiducial region, /;-[ > 0.1 GeV, 
and for full acceptance, pj > 0 GeV. The results are presented in //-intervals of 0.1 unit width. Due to the 
excellent //-resolution of the tracklets, as seen from Fig. 4, migration of tracklets between neighbouring 
bins is negligible. 
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The fully corrected, per-event charged-particle pseudorapidity distributions arc calculated according to 

dNch _ 1 E AMr(O' Zvtx, >l)C(Q, Zv tx , ??) , 5 ) 

d ?7 A rj Yj A^evtCZvtx) 

where A/V tr indicates either the number of reconstructed pixel tracks or two-point tracklets, M ev t(Zvtx) is 
the number of analysed events in the intervals of the primary vertex along the z direction, and the sum 
in Eq. (5) runs over primary vertex intervals. The number of primary vertex intervals varies from seven 
for |//| < 2.2 for two-point tracklets and |//| < 2 for pixel tracks to two at the edges of the measured 
pseudorapidity range of |//| < 2.7 for two-point tracklets and \r/\ < 2.5 for pixel tracks respectively. The 
primary vertex intervals used in the analysis are chosen such that C{0,z vix ,rf) changes by less than 20% 
between any pair of adjacent z vtx intervals. 

Figure 6 shows the effect of the applied correction for all three methods. The left panels shows the MC 
simulation results based on Hijing. The distribution of generated primary charged particles is shown by 
a solid line and the distributions of reconstructed tracks and tracklets are indicated by markers in the 
upper left panel. The lower left panel shows the ratio of reconstructed distributions to the generated 
distribution. Among the three methods, the corrections for Method 1 are the smallest, while the pixel 
track method requires the largest corrections. The structure of the measured distribution for the pixel 
track method around 77 = ±2 is related to the transition between the barrel and endcap regions of the 
detector. The open markers in the right panel of Fig. 6 show the reconstructed distribution from the data 
and the tilled markers are the corresponding distribution for the three methods after applying corrections 
derived from the simulation. The lower panel shows the ratio of the results obtained from Method 2 and 
the pixel track method to that obtained using Method 1. The three methods agree within 2% in the barrel 
region of the detector and within 3% in the endcap region. This agreement demonstrates that the rejection 
of fake track or tracklets and the correction procedure are well understood. For this paper. Method 1 is 
chosen as the default result for dA^h/dr/, Method 2 is used when evaluating systematic uncertainties, and 
the pixel track method is used primarily as a consistency test, as discussed in detail below. 


7 Systematic uncertainties 


The systematic uncertainties on the d/V ch /d?/ measurement arise from three main sources: inaccuracies in 
the simulated detector geometry, sensitivity to selection criteria used in the analysis including the residual 
contributions of fake tracklets and secondary particles, and differences between the generated particles 
used in the simulation and the data. To determine the systematic uncertainties, the analysis is repeated 
in full for different variations of parameters or methods and the results are compared to the standard 
Method 1 results. A summary of the results are presented in Table 1. 

The uncertainty due to the simulated detector geometry arises primarily from the details of the pixel 
detector acceptance and efficiency. The locations of the inactive pixel modules are matched between the 
data and simulation. Areas smaller than a single module that are found to have intermittent inefficiencies 
are estimated to contribute less than 1.7% uncertainty to the final result. This uncertainty has no centrality 
dependence, and is approximately independent of pseudorapidity. 

The amount of inactive detector material in the tracking system is known with a precision of 5% in the 
central region and up to 15% in the forward region. In order to study the effect on the tracking efficiency, 
samples generated with increased material budget are used. The net effect on the final result is found to 
be 0.5-3% independent of centrality. 
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Uncertainties due to tracklet selection cuts are evaluated by independently varying the cuts on |Ap| and 
|A0| up and down by 40%. The effect of these variations is less than 1%, except at large values of \p\ 
where it is 1.5%, and has only a weak centrality dependence. 

The systematic uncertainty due to applying the pj reweighting procedure to the generated particles is 
taken from the difference in dN c h/dr/ between applying and not applying the reweighting procedure. The 
uncertainty is less than 0.5% for |//| <1.5 and grows to 3.0% towards the edges of the // acceptance. The 
uncertainty has a centrality dependence because the pj distributions in central and peripheral collisions 
are different. 

Tracklets are reconstructed using Method 1 for particles with pj > 0.1 GeV. The unmeasured region of 
the spectrum contributes approximately 6% to the final d/V t h/d?/ distribution. The systematic uncertainty 
on the number of particles with pj <0.1 GeV is partially included in the variation of the tracklet A (p se¬ 
lection criteria. An additional uncertainty is evaluated by varying the shape of the spectra below 0.1 GeV. 
This uncertainty is estimated to be as much as 2.5% at large values of \rj\ and has a weak centrality 
dependence. 

To test the sensitivity to the particle composition in Hi.iing, the fraction of pions, kaons and protons in 
Hijing are varied within a range based on measured differences in particle composition between pp and 
Pb+Pb collisions [60, 61]. The resulting changes in dV c h/d ?7 are found to be less than 1% for all centrality 
intervals. 

Systematic uncertainties due to the fake tracklets are estimated by comparing the results of the two tracklet 
methods. The differences in the most central collisions are found to vary with pseudorapidity from 1.5% 
in the barrel region to about 2.5% at the ends of the measured pseudorapidity range. 


Source 

60-90% centrality 
barrel endcap 

0 -1% centrality 
barrel endcap 

Inactive modules 


1.7% 

1 . 

7% 

Extra material 

0.5% 

3.0% 

0.5% 

3.0% 

Tracklet selection 

0.5% 

1.5% 

0.5% 

1.5% 

pj reweighting 

0.5% 

0.5% 

0.5% 

3.0% 

Particles with pj < 0.1 GeV 

1 .0% 

2.5% 

1 .0% 

2 .0% 

Particle composition 


1 .0% 

1 . 

0 % 

Contribution of fake tracklets 

1.5% 

2 .0% 

1.5% 

2.5% 

Event selection efficiency 

5.0% 

6 .0% 

0.5% 

0.5% 

Total 

5.7% 

7.9% 

2.9% 

5.9% 


Table 1: Summary of the various sources of systematic uncertainty and their estimated impact on the dV c h/d;/ 
measurement in central (0-1%) and peripheral (60-90%) p +Pb collisions. 


The uncertainty associated with the event selection efficiency for the fiducial class of p +Pb events is eval¬ 
uated by defining new cen fi' a lity ranges after accounting for an increase (decrease) in the efficiency 
by 2% and repeating the full analysis. This resulting change of the dV C h /d// distribution is less than 0.5% 
in central collisions; it increases to 6% in peripheral collisions. 

The uncertainties from each source were evaluated separately in each centrality and pseudorapidity to 
allow for their partial or complete cancellation in the ratios of d/V c h /dp distributions. The impact in 
different regions of pseudorapidity and centrality are shown in different columns of Table 1 . Uncertainties 
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coming from different sources and listed in the same column are treated as independent. The resulting 
total systematic uncertainty shown in the lower line of the table is the sum in quadrature of the individual 
contributions. 


8 Results 


Figure 7 presents the charged-particle pseudorapidity distribution dA^h/d?/ for /;+Pb collisions at y.v KN . 
=5.02 TeV in the pseudorapidity interval |//| < 2.7 for several centrality intervals. The left panel shows 
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Figure 7: Charged-particle pseudorapidity distribution measured in several centrality intervals. Left: diV C h/dp for 
charged particles with p-\ > 0.1 GeV. Right: d/V ch /<P/ for charged particles with p-\ > 0 GeV. Statistical uncer¬ 
tainties, shown with vertical bars, are typically smaller than the marker size. Shaded bands indicate systematic 
uncertainties on the measurements. 


the dfVch/dp distribution measured in the fiducial acceptance of the ATLAS detector, detecting particles 
with pj > 0.1 GeV. The results for the d/V c h/d?/ distribution with pj > 0 GeV are shown in the right 
panel of Fig. 7. The charged-particle pseudorapidity distribution increases by typically 5%, consistent 
with extrapolation of spectra measured in pp collisions to zero pt [40]. At the edges of the measured 
pseudorapidity interval, it increases dN c h/dp by 11%. 

In the most peripheral collisions with a centrality of 60-90%, the d/V ch /dp distribution has a doubly- 
peaked shape similar to that seen in pp collisions [40, 62]. In collisions that are more central, the shape 
of d/V ch /d?/ becomes progressively more asymmetric, with more particles produced in the Pb-going dir¬ 
ection than in the proton-going direction. To investigate further the centrality evolution, the dNch/dp 
distributions in each centrality interval are divided by the dA c h/dp distribution for the 60-90% interval. 
The results are shown in Fig. 8, where the double-peak structure disappears in the ratios. The ratios are 
observed to grow nearly linearly with decreasing pseudorapidity, with a slope whose magnitude increases 
from peripheral to central collisions. In the 0-1% centrality interval, the ratio changes by almost a factor 
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Figure 8: Ratios of dA^h/d;/ distributions measured in different centrality intervals to that in the peripheral (60-90%) 
centrality interval. Lines show the results of second-order polynomial fits to the data points. 


of two over the measured p-range. The greatest increase in multiplicity between adjacent centrality inter¬ 
vals occurs between the 1-5% and 0-1% intervals. Averaged over the //-interval of the measurement, the 
d/V c h/d// distribution increases by more than 25% between the 1-5% and 0-1% intervals. 
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Figure 9: Ratios of d/V c |,/d// obtained using alternative centrality definitions to the nominal results presented in this 
paper as a function of r\ for the 0-1% and 60-90% centrality bins. 


In addition to the results presented in Figs. 7 and 8, the dN c h/drj measurement is repeated using the altern¬ 
ative definitions of the event centrality variables defined in Sect. 5. Figure 9 demonstrates the sensitivity 
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of the measured dA c n/d// to the choice of centrality variable by showing the ratios of the dA^h/d;/ distri¬ 
butions in the most central and most peripheral intervals under the + and 2£'^ 6< ' ,?cm ' <44 centrality 

definitions to those obtained with the nominal definition. Using the XE-,- 4 centrality definition, 
the diV c h/ dr] distributions change in an approximately ^-independent fashion by -3% and +3% for the 0- 
1% and 60-90% intervals, respectively. The d/V ch /dt/ distributions in the other centrality intervals change 
in a manner that effectively interpolates between these extremes. As a result, the increase in dN c h/dr/ 
between the most peripheral and most central collisions would be reduced by 6 % relative to the nom¬ 
inal measurement. Using the symmetric, cen t ra ]ity definition, the d/V c h/d?/ distribution in 

each interval changes in an 77 -dependent way such that the ratio is consistent with a linear function of 77 . 
The change is at most 6 % at the ends of the 77 range in the most central and most peripheral centrality 
intervals, and smaller elsewhere. Thus, for the symmetric centrality selection the ratios in Fig. 8 for the 
0-1% bin would increase by 9% at 77 = 2.7, and decrease by 6 % at 77 = -2.7. Generally, the alternative 
centrality definitions considered in this analysis yield no qualitative and only modest quantitative changes 
in the centrality dependence of the d/V t h/d?/ distributions. These variations should not be considered a 
systematic uncertainty on the dN c h/dij measurement but do indicate that the particular centrality method 
used in the analysis must be accounted for when interpreting the results of the measurement. 

Figure 10 shows a comparison, where possible, of the measurements presented in this paper to results 
from the ALICE experiment [31] using a centrality definition that is based on the detector covering the 
pseudorapidity region -5.1 < 77 < -2.8, similar to the ^E^-bascd selection used in this measurement. 
The ATLAS results for 0-1% and 1-5% centrality intervals are combined to match the ALICE experiment 
result for 0-5% interval. Similarly, the 20-30% and 30—40% intervals are combined to match the ALICE 
experiment result for 20-40% interval. The results from the two experiments are consistent with each 
other. 
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Figure 10: Charged-particle pseudorapidity distribution dAAh/tl// measured in different centrality intervals compared 
to similar results from the ALICE experiment [31] using the “VBA” centrality selection. The ATLAS centrality 
intervals have been combined, where possible, to match the ALICE centrality selections. 


1 1 1 I 1 1 1 

ATLAS 

p+Pb, 1 pb 
\ s NN = 5.02 TeV 
= 0.465 


0-5% 


ATLAS 

ALICE 


= u.4bo 

5-10% •**•... 

..-TTTt... ''• Tt *?TTTTTT“TTt T ..,. 

10 - 20 % 

20-40% ’’’’Wftmwmtwm,,,.]’***. 

«... 

40 60°/ 


17 




9 Particle multiplicities per participant pair 


A common way of representing the centrality dependence of particle yields in A+A and p+A collisions 
is by showing the yield per participant or per participant pair, (A part )/2, which is determined for each 
centrality interval and each geometrical model as shown in Fig. 3. Figure 11 shows dA c h/d?/ per parti¬ 
cipant pair for the most central and most peripheral intervals of centrality measured as a function of ij for 
three different models of the collisions geometry: the standard Glauber model and the GGCF model with 
oj, t = 0.11 and 0.2 in the top, middle and lower panels, respectively. The results for the most peripheral 



0 

Figure 11: Charged-particle pseudorapidity distribution d/V c h/cl// per pair of participants as a function of;/ for 0-1% 
and 60-90% centrality intervals for the three models used to calculate A part . The standard Glauber calculation is 
shown in the top panel, the GGCF model with = 0.11 in the middle and o> <T = 0.2 in the lowest panel. The bands 
shown with thin lines represent the systematic uncertainty of the d/V L |,/d// measurement, the shaded bands indicate 
the total systematic uncertainty including the uncertainty on (Apart). Statistical uncertainties, shown with vertical 
bars, are typically smaller than the marker size. 


(60-90%) centrality interval, shown with circles, are similar between all three panels. This is due to relat¬ 
ively small difference between the calculations of (Apart) for Glauber and GGCF models in this centrality 
interval. The shape of the distribution indicates more abundant particle production in the proton-going 
direction in comparison to the Pb-going. This can be explained by the higher energy of the proton com¬ 
pared to the energy of a single nucleon in the lead nucleus in the laboratory system. In the most central 
collisions (0-1%), shown with diamond markers in all three panels, this trend is reversed. Conversely, 
the magnitude of dA c h/d// per participant pair strongly depends on the geometric model used to calculate 
(Apart). The point at which the central and peripheral scaled distributions cross each other also depends 
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on the choice of geometric model. 


Figure 12 shows the dA^h/dr? distribution per participant pair as a function of (N part ) for the three differ¬ 
ent models of the collisions geometry. Since the charged-particle yields have significant pseudorapidity 
dependence, dA r c h/d?7/((A r part)/2) is presented in five rj intervals including the full pseudorapidity interval, 
-2.7 < ij < 2.7. In the region 0 < ij < 1, the d/V ch /d// distribution is consistent with an empirical fit to 
inelastic pp data that suggest dA^h/d// increases with centre-of-mass energy, -\[s, as (or .y 010 ^ [16]. 



Figure 12: Charged-particle pseudorapidity distribution d;V L . h /d;/ per pair of participants as a function of (A pal1 ) in 
several ^-regions for the three models of the geometry: the standard Glauber model (top panel), the GGCF model 
with Wq. = 0.11 (middle panel) and GGCF with oj it = 0.2 (bottom panel). The open boxes represent the systematic 
uncertainty of the d;V c i,/d// measurement only, and the width of the box is chosen for better visibility (they are not 
shown for -1.0 < rj < 0 and 0 < rj < 1). The shaded boxes represent the total uncertainty (they are shown only 
on -2.7 < rj < 2.7 interval for visibility) which is dominated by the uncertainty of the (N p - dn ) given in Table 4 and 
Fig. 3. This uncertainty is asymmetric due to the asymmetric uncertainties on (N pan ). The statistical uncertainties 
are smaller than the marker size for all points. 


The d/Veh / d ///((/ V p ait )/ 2 ) values from the standard Glauber model are approximately constant up to (N pan ) ~ 
10 and then increase for larger (N part ). This trend is absent in the GGCF model with oj <t = 0 . 11 , which 
shows a relatively constant behaviour for the integrated yield divided by the number of participant pairs. 
The d / V t h /dr//(( A^pait )/ 2 ) values from the GGCF model with to a = 0.2 show a slight decrease with ( A ^ par t ) 
in all rj intervals. 

The presence or absence of (A^art) scaling does not suggest a preference for one or another of the geo¬ 
metric models. Flowever, this study emphasises that considering fluctuations of the nucleon-nucleon 
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cross-section in the GGCF model may lead to significant changes in the V part scaling behaviour of the 
p+Pb d/V ch /d// data and, thus, their interpretations. 


10 Conclusions 

This paper presents a measurement of the centrality dependence of the charged-particle pseudorapidity 
distribution, d/V c h /dp, measured in approximately 1 pfr 1 of p+Pb collisions at a nucleon-nucleon centre - 
of-mass energy of y.v N . N . =5.02TeV collected by the ATLAS detector at the LHC. The fully coiTected 
measurements arc presented for the fiducial acceptance of the ATLAS detector ( p-y >0.1 GeV) and in the 
full acceptance (pj > 0 GeV). The d/V c h/dp distributions arc presented as a function of pseudorapidity 
over the range -2.7 < p < 2.7 and as a function of collision centrality for the 0-90% p+Pb collisions. The 
centrality is characterised using the energy deposited in the forward calorimeter covering -4.9 < p < -3.1 
in the Pb-going direction. 

The shape of d/V c h /dp evolves gradually with centrality from an approximately symmetric shape in the 
most peripheral collisions to a highly asymmetric distribution in the most central collisions. The ratios of 
dVch/dp measured in different centrality intervals to the dfV c h/dp distribution in the most peripheral inter¬ 
val arc approximately linear in p with a slope that is strongly dependent on centrality. It is noteworthy that 
the greatest increase in charged-particle multiplicity between successive centrality bins occurs between 
the 1-5% and 0-1% centrality bins. 

The results arc also interpreted using models of the underlying collision geometry. The average number 
of participants in each centrality interval, (V part ), is estimated using a standard Glauber model Monte 
Carlo simulation with a fixed nucleon-nucleon cross-section, as well as with two Glauber-Gribov colour 
fluctuation models which allow the nucleon-nucleon cross-section to fluctuate event-by-event. The V part 
dependence of d V t h / d p / (((V p a rt ) / 2) is found to be sensitive to the modelling of the p+Pb collision geo¬ 
metry, especially in the most central collisions: while the standard Glauber modelling leads to a strong 
increase in the multiplicity per participant pair for collisions in the centrality range (0-30)% the GGCF 
model produces a much milder centrality dependence. 

These results point to the importance of understanding not just the initial state of the nuclear - wave func¬ 
tion, but also the fluctuating nature of nucleon-nucleon collisions themselves. 
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Appendix: Glauber model analysis 


The PHOBOS Glauber MC program [53] is used to perform the standard Glauber model calculations 
used in this analysis. The Pb nucleon density is taken to be a Woods-Saxon distribution with radius 
and skin depth parameters, R = 6.62 fm and a = 0.546 fm [63], respectively. The nucleon-nucleon 
inelastic cross-section is taken to be 70 mb. The resulting probability distribution, P(N part ), of the number 
of participating nucleons N pan - nucleons that undergo at least one hadronic scattering during the p +Pb 
collision - is shown in Fig. 13. 



Figure 13: Glauber and GGCF N part distributions for 5.02 TeV p+Pb collisions obtained from one million simulated 
events each. The inset shows the GGCF PhOt^) distributions for oj :r = 0.11 and 0.2. 


The GGCF model is implemented in a modified version of the PHOBOS MC program. Following 
Ref. [36], the probability distribution to find the nucleons in a configuration having a nucleon-nucleon 
scattering cross-section, cr, is taken to be 


P{cr) = p 


cr 

cr + cr 0 


exp 


(cr/cro - l) 2 } 

n 2 J' 


( 6 ) 


Here, p is a normalisation constant, Q controls the width of the P(<x) distribution, and cr 0 determines the 
configuration-averaged total cross-section cr tot = (cr). The inelastic cross-section, cr NN ,, is taken to be a 
constant fraction, A, of the total cross-section [37] so the probability distribution of <r SK is given by 

^h(ct nn ) = - { P(cr m /A). (7) 


The values used in this analysis for Q, cro, cr tot and A corresponding to oj it = 0.11,0.2 are shown in 
Table 2. The first, earlier analysis yielding (o a = 0.11 [36] assumed cr tot = 86 mb, consistent with the 
Donnachie and Landshoff [64] parameterisation of cr tot ( s ). The second analysis yielding o)< r = 0.2 used 
an updated measurement of the pp total cross-section at the LHC [65] to set cr tot = 94.8 mb. However, 
modifying the parameters for the oj a - 0.11 case to be consistent with this improved knowledge of cr tot 
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produces a negligible change in the resulting P(cr) distribution. The values for A are chosen to produce 
the above-quoted nucleon-nucleon inelastic cross-section of 70 mb. The GGCF distributions 

are shown in the inset of Fig. 13, while the resulting P(A parL ) distributions are shown in the main panel of 
the figure. 

To connect an experimental measurement of collision centrality such as X F [ p b to the results of the Glauber 
or GGCF Monte Carlo simulation, a model for the A part dependence of the X-Fx b distribution is required. 
The usual basis for models previously applied to A+A and p/d+A collisions is the WN model [12], which 
predicts that the average increases proportionally to /V parl with the proportionality constant equal to 
one half the corresponding average FCal X^r in pp collisions. 

Under the WN model, the X E^ b distribution for fixed Apart would be obtained from a A parL -fold convo¬ 
lution of the corresponding distribution in pp collisions. This convolution is straightforward if the X^t 
distribution in pp collisions is described by a gamma distribution [66] 

gamma (Z£t;M) = (8) 

since gamma distributions have the property that an /V-fold convolution of a gamma distribution with 
parameters k and 6 yields another gamma distribution with the same 6 and a modified k parameter, k! = 
Nk. 

Attempts to tit the measured Y,E bb distribution using pure WN-convolved gamma distributions and the 
Glauber /V part distribution yield unphysical results for the nucleon-nucleon parameters, ko and 9o, when 
those parameters are free parameters of the fit. In particular, ko is less than unity, which implies a X-Ft 
distribution that increases with decreasing X^r faster than e anc j q {) j s unrealistically large. The 

resulting nucleon-nucleon X^r distribution is also inconsistent with that measured in pp collisions [67]. 
The poor behaviour of the WN model is primarily due to the difference in shape between the Glauber 
Apart distribution and the measured XF- p h distribution. To improve the description of the measured X E? b 
distribution, a generalisation of the WN model is implemented that parameterises the A part dependence of 
the k and 9 parameters of the gamma distribution as 

k (A part ) - ko + k\ (A pan — 2), 

e{N VArl ) = 6 {) + 6 ] log(A part - l). (9) 

For k\ = & 0/2 and 9\ = 0, this model reduces to the WN model. The log (A part - I j term allows for a 
possible variation in the effective acceptance of the FCal due to an Ap art -dependent backward shift in the 
p +Pb centre-of-mass system [68]. 

To limit the number of free parameters when fitting the ^Ej b distribution, ko and 9o are obtained by 
fitting the detector-level X A p distributions in Pythia6 and Pythia8 pp simulations. These simulations 


Parameter 

(x>(j — 0.11 

(Do- = 0.2 

Q 

0.55 

1.01 

cro [mb] 

78.6 

72.6 

cr tot [mb] 

86 

94.8 

A 

0.82 

0.74 


Table 2: Parameters used in the parameterisation of the GGCF P(cr tot ) distribution. 
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have been shown to give a reasonable description of the corresponding pp collision data at \fs = 7 TeV 
[67] although they both slightly under-predict the average forward transverse energy. The contribution 
of electronic noise to the simulated distribution was determined by examining the distribution in 

empty beam bunch crossings in data. In the fit, the gamma distributions were convolved with the effects 
of this noise before comparison with the data. The Pythia8 fit results, ko = 1.40 and 0o = 3.41, are used 
for the default analysis. The Pythia6 fit results, ko = 1.23 and 6q = 2.68, are used to evaluate systematic 
uncertainties. 

The measured X E^ 3 distribution is fitted with a distribution produced by summing the fVp art -dependent 
gamma distributions, after weighting them by P(/V pa , t ) and including an additional convolution to account 
for electronic noise. The model distribution is also re-weighted to properly describe the ££;!j! b -dependent 
event selection efficiency in the data, which is estimated using the Pythia MC samples under the assump¬ 
tion that the p +Pb inefficiency for a given ' s same as that in pp collisions. Results are shown 

in Fig. 14 for the Glauber model and the two GGCF models. The fits provide a good description of the 



Figure 14: Top panel: Measured distribution compared to Glauber (solid), GGCF with oj lr = 0.2 (long 

dashed), and GGCF with o>, r = 0.11 (short dashed) fits. Lower panels: ratios of Glauber and GGCF fit distributions 
to the data distribution. 

distribution for 2^x b <100 GeV for all three geometric models, although at higher the Glau¬ 
ber fit describes the data better. The deviations of the GGCF fits from the data become significant near 
= 120 GeV; the fraction of the total X E ] ’ h distribution above this value is approximately 0.1%. 

The parameters k\ and 0\ are obtained from fixing ko and do and fitting the distribution to the data. 
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They are presented in Table 3 along with the ratios ki/ko for each of the geometric models. Pure WN 
behaviour would correspond to k\ /ko = 0.5 and 6 1 = 0. The results indicate substantial deviations from 
WN behaviour for the Glauber and GGCF coa- = 0.2 fits, while the GGCF oj ir = 0.11 fit yields both a 
k\/ko that is close to 0.5 and a small 0\ . The success of the above-described fitting procedure in describing 
the measured J]Ej b distributions using parameterisations of the /V parl dependence of the J]Ej b response 
with only two free parameters is due to the similarity of the shapes of the P(N parl ) distribution and the 
measured X-Fj 5 distributions. However, the pronounced knee in the Glauber P(N par t ) distribution requires 
more non-linearity in the N parl dependence of the XF p h response. In contrast, the lack of such a feature 
in the GGCF P(N pa n) distributions allows a simpler description of the measured X F- p h distribution. 

The results of the fit procedure described above provide a data-driven estimate of the total p+ Pb event- 
selection efficiency. For the default PYTHiA8-based results, the integral of the simulated XF p b distribution 
is 2% higher than that in the data. The deficit in the data is concentrated at low XFj h values, consistent 
with losses due to event selection. However, a detailed analysis of the residual differences between the 
best fit and measured distributions indicates an excess of very low X-Fj 13 events i n the data, which varies 
from ~ 0% for the Glauber fit to 1.8% for the GGCF tit with oj ir = 0.2. These may arise from residual 
diffractive or photo-nuclear collisions and are considered background. For the purpose of defining X F. p h 
centrality intervals, this background effectively increases the event-selection efficiency by adding events 
that are all in the 90-100% centrality interval. For Pythia6, the X Ej fits using the default model yield 
a total efficiency of 97% and up to 1% background. The alternative models for k(jtyp art ) and d(/V parl ) 
yield a similar total efficiency of 98% and a background rate as high as 2%, a rate that is compatible with 
independent estimates of the rate for collisions involving diffractive excitation of the proton to pass the 
applied event selections. Based on these results, the total effective efficiency including background is then 
taken to be 98% and the uncertainty is conservatively estimated to be 2%. 

The (./Vp ar t> values are obtained for each of the centrality intervals using the results of the fits to the X F. p h 
distributions. The ( N part ) results along with the total systematic uncertainties, which are described below, 
are shown in Fig. 3 and listed in Table 4. 

To obtain systematic errors on N paTt in each centrality interval, the maximum positive and negative frac¬ 
tional variation in (N pan ) away from the default results is determined for different classes of variations, 
detailed below. 

To evaluate the impact of the total event selection uncertainty, new centrality intervals are chosen assum¬ 
ing a total efficiency of 100% and 96% and the complete analysis is repeated. To account for possible 
inaccuracies in the PYTHiA8-simulated dEj/dij in the region of the FCal acceptance, the analysis is re¬ 
peated separately under ±10% re-scalings of Pythia8 XF^ 13 va ' ucs commensurate with the scale of the 
data-PYTHiA8 differences observed in Ref. [67], Other variations for which the complete analysis is re¬ 
peated are: (i) using the Pythia6 event generator to fix ko and Go, (ii) alternative models for k (/V pail ) and 

6 (Npan), (iii) ±5 mb changes in <x N . N ., and (iv) variations in the parameters of the nuclear density distri¬ 
bution. For the model uncertainty, two alternative parameterisations for k (/V part j and f) (/V pa n) are used. 
One of these kept 6 constant, 0 (/V pait ) = Go while allowing for a quadratic dependence of k on /V part . The 

other included both a quadratic term in k^N par t ) and the logarithmic term in d(<V part j but fixed k\ = ko/2 
to reduce the number of free parameters. 

The resulting maximal variations are then summed in quadrature over the different classes separately 
for the positive and negative variations to produce the uncertainties listed in Table 4. The uncertainties 
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Glauber model 

ki 

ki/ko 

9i 

Standard Glauber 

0.425(2) 

0.304(1) 

+ 1.32(1) 

GGCF ov = 0.11 

0.901(3) 

0.643(2) 

+0.074(4) 

GGCF air = 0.2 

1.139(3) 

0.813(2) 

-0.209(2) 


Table 3: Optimal fit parameters and uncertainties obtained from fits to the measured ££’. l | >b distribution. Uncertain¬ 
ties on the fit parameters are shown in parenthesis. 


Centrality 

Glauber 

GGCF oj tT 

= 0.11 

GGCF uv = 0.2 

60-90% 


'+5%) 

-8 %) 


+5%\ 

-5%) 

3. 4 tl( 

+8%) 

- 5 %; 

40-60% 

7 +S:d 

4 - 6 % \ 
-8%) 

6-6^ 

+6%\ 

-6%) 


+8%) 

- 5 %; 

30^10% 

9 8 +0 - 6 
*-0.6 

4 - 6 % \ 
-6%) 

9.2«j( 

+6%\ 

-6%) 

8 . 9 ;“ ( 

+ 7 %) 

- 5 %; 

20-30% 

i i 4+0.6 
11 -0.6 

(+6%\ 

\-6%) 

11 - 2^7 

( +6%\ 
V-6%) 


(+6%) 

1-6%; 

10-20% 

13.0+Q7 

(+6%\ 

\-6%) 

1 4 7+0.8 
J -0.8 

(+6%) 

\-7%/ 

14 l +0 - 9 
-0.8 

(+6%\ 

1-6%; 

5-10% 

14 - 6 -a8 

( +8%\ 
\-6%) 

16 5 + 1 ° 

ru.J i o 

(+6%\ 

V-6%) 

17 4 +hl 

-l.t 

/ + 7 %) 
1-6%; 

1-5% 


f+11%) 

1-6% ) 

19 5 +1 - 3 
i *‘ -1.3 

(+ 7 %) 

1 - 7 %; 

21 - 4 -2 : 4 o 

(+i%\ 

1-9%; 

0-1% 

"U-fl 

f+15%) 

1-5% ) 


(+1%\ 

i-8%; 

27 - 4 -4:o< 

%6% \ 
,-ie%; 

0-90% 

8 - 4 td 

4 - 6 % \ 
k-5%) 

8 -5iS:4 

+6%) 

- 5 %; 

8.6^5 ( 

+6%) 

- 5 %; 


Table 4: (/V p;lrl ) values for centrality intervals used in this analysis together with asymmetric systematic uncertainties 
shown as absolute as well as relative uncertainties. 


for most centrality intervals are dominated by comparable contributions from the choice of model para- 
meterisation, differences between the Pythia6 and Pythia8 fits, and uncertainties in cr NN . For the 
40-60% and 60-90% interval the uncertainty in event selection efficiency has a contribution to the _/V part 
systematic uncertainty that is similar in magnitude to these other three sources. 
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